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tional  to  the  shock  velocity  in  the  explosive.  For  true  one-dimensional  situations, 

Che  area  dependant  power  requirement  is  indistinguishable  from  the  constant  energy  (per 
unit  area)  criteria. 

An  Important  spinoff  from  this  investigation  is  Che  forsulaCion  of  a  scheme  to  predict 
Che  minioum  shock  pressure  required  to  detonate  an  explosive.  This  scheme  is  based  on 
E.  R.  Flttgerald's  crystal  lattice  fracture/disintegration  criteria.  Some  pertinent 
results  are  given  in  this  report,  but  Che  concept  and  analysis  are  being  documented  in. 
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I.  INTRODUCTION 


Criteria  for  high  explosive  detonation  via  projectile  impact  are  important 
from  both  safety  and  lethality  considerations.  The  study  of  both  impact  and 
explosive  detonation  encompasses  some  of  the  most  complex  phenomena  encoun~ 
tered  in  the  physical  sciences.  Consequently,  each  area  (high  velocity  Impact 
and  explosive  detonation)  has  been  the  subject  of  numerous  conferences, 
symposiums,  and  technical  papers. 

Since  projectile  Induced  detonation  criteria  has  such  important  lethality, 
safety,  and  survivability  implications,  a  technical  literature  information 
search  was  initiated.  K  primary  objective  was  to  acquire  insight  into  the 
physical  phenomena  causing  detonation  via  impact.  This,  in  turn,  could  lead 
to  better  utilization,  or  even  enhancement,  of  present  empirical  predictive 
schemes  employed  in  various  lethality  codes.  Reference  1  employs  part  of 
this  predictive  methodology  (the  Jacobs-'Roslund  equation)  in  the  analysis  of 
experimental  data. 

This  information  search  has  yielded  five  sets  of  experimental  data  which 
have  been  reexamined  in  this  investigation.  These  data  have  been  derived  by 
subjecting  bare  unprotected  explosive  samples  to  normal  impact  of  flat-faced 
projectiles  as  schematically  shown  in  Figure  1.  The  basic  information,  which 
came  from  diverse  sources,  is  plotted  in  Figure  2. 

The  information  search  is  cootinuiug  and  appropriate  data  are  being  reex¬ 
amined.  However,  enough  analysis  has  been  done  to  indicate  that  enhancements 
or  supplements  to  certain  portions  of  the  current  empirical  detonation  predic¬ 
tion  methodology  may  be  feasible  or  even  necessary.  In  any  event,  the  current 
methodology  (for  prediction  of  explosive  detonation  via  projectile  velocity 
and  size)  could  be  checked  by  incorporating  the  "vcritical  area’  (A^r)  concept 
suggested  by  Henry  Moulard  [2]  and  [3]. 

Specifically,  during  the  present  investigation,  some  empirical  relation¬ 
ships  were  found  which  are  simple  functions  of  known  impact  shock  variables 
and  either  or  A,  the  projectile  cross-section  arcs. 

These  relations,  involving  cross-section  area  dimensions,  can  be  useful  in 
projectile  size/detonation  prediction  criteria,  for  lechallty/survivability 

assessment  codes. 

Even  though  the  physical  meaning  of  certain  empirical  relations  is  not 
crystal  clear,  their  uniform  and  consistent  data  correlation  suggests  that 
they  are  natural  parameters  associated  with  fundamental  physical  phenomena. 
While  searching  for  some  explanation  of  the  strong  critical  area  effect,  its 
connection  with  an  energy  rate  or  power  input  was  found. 

The  present  investigation  strongly  indicates  that  a  minimum  (or  critical) 
area  dependent  shock  energy  input  rate  (or  power)  is  necessary  for  detonation 
initiation.  This  is  analogous  to  results  reported  by  investigators  of  gaseous 
explosive  detonation.  This  mlmumum  amount  of  power  input.  Joules/microsecond 
(J/usec),  (Dust  be  sustained  for  a  minimum  amount  of  time  (ter) »  which  is  de¬ 
pendent  on  the  overall  energy  input  requirements. 


It  is  ahovra  that  the  area  dependent  energy  rate  la  Inversely  proportional 
to  the  shock  velocity  in  the  explosive*  An  obvious  question  is:  "What  is  the 
ttinimum  shock  velocity  required  to  initiate  reactions  leading  to  detonation?" 
This  question  led  to  the  results  described  in  Appendix  G.  In  this  appendix,  a 
crystal  lattice  fracture  or  disintegration  criteria  is  applied  to  compute  the 
mlnioum  shock  velocity  and  pressure  required  to  initiate  detonation*  More 
work  remains  to  be  done  with  respect  to  this  concept,  since  it  appears  to  be 
fundamentally  significant* 

Tima  visa  integration  of  the  area  dependent  energy  input  rate  expression 
produced  en  expression  for  the  total  high  shock  energy  input  as  a  function  of 
time*  This  expression  is  employed  to  demonstrate  that  when  the  projectile  and 
explosive  cross  section  areas  are  sufficiently  larger  than  Chen  the  well 

known  (area  independent)  constant  energy  per  unit  area  detonation  criteria 
will  logically  appear  as  the  dominant  parameter*  This  situation  is  prevalent 
in  the  one^lmcnslonal  tests  and  it  may  obscure  the  important  energy  rate  (or 
power)  input  role* 

The  seven  appendices  of  this  report  provide  certain  background  information 
with  respect  to: 

A,  B  -  Shock  Phenomenon  and  Data 

C  -  The  Critical  Energy  Criteria 

D  -  The  Critical  Area  Criteria 

E  -  Current  Projectile  Impact  Prediction  Criteria 

F  '■  Explosive  Rod  Diameter  Effect 

G  -  Crystal  Lattice  Fracture  and  Detonation  Criteria 
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II.  ANALYSIS 


A.  Critical  Area  Acr  Relations 

Moulard  (2)  clearly  demonstrated  hla  point  with  respect  to  A^-r  for  an 
explosive  which  conformed  to  the  constant  energy  criteria  (Pj  'Jpgjj  ^cr  “ 
Const).  The  constant  energy  relation  [4]  was  developed  from  an  examination  of 
one  dimensional  thin  metallic  foil  planar  Impact  detonation  data  reported  In 
[5].  Moulard  points  out  that  small  projectile  or  fragment  Impact  la  a  three- 
dimensional  phenomenon,  even  If  the  contact  surfaces  are  planar.  Moulard 
showed  that  a  minimum  shock  loaded  Interface  contact  area  (Agr)  Is  Important 
In  addition  to  the  well  known  shock  pressure  (Pg)  and  duration  (t^r)  condi¬ 
tions  from  the  constant  energy  per  unit  area  relation: 


Pq  Un  c  "Cl  -  Constant  (1) 

®  Pex  '■cr  ^ 

^Pex  “  Particle  velocity  of  the  explosive  behind  the  shock  front. 

A  dimensional  analysis  of  Eq.  (1)  reveals  that: 


P.  Up.,  'cr  -  ^ 

L  m  ^O’^ce 

"  L  ”  Length 


Work 

Area 


(2) 


This  shows  that  force  per  unit  length  could  also  be  important.  Was 
there  some  other  combination  of  variables  having  dimensions-  of  force  per  unit 
length  which  was  constant  also?  Subsequent  work  with  Moulard 's  [2]  data  for 
the  Comp-B3  explosive  (Table  1  and  Appendix  D)  revealed  empirically  that: 


■  C2  •  Constant 

_  /f  \l  _  Work 

”  \L2y  “  Area 

_  /fN  _  Force 

yL  J  Length 


f3) 


Note  Chat  Moulard  [3]  presented  correlation  logarithmic  plots  of 
Pg  versus  A^f.  Equation  (3)  explicitly  delineates  a  plausible  Pg  and  Agf 
relationship. 


Division  of  Eq.  (1)  by  Eq.  (3)  yields: 
^Pex  ^cr  _  ^1  _  r. 

■  cl  ■  '=3 


■  Dimensionless  (4) 

Constant 

The  [2]  data  also  confirmed  Eq.  (4). 
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Would  any  of  Che  ocher  available  small  projecClle  Impacc  detonacion 
daCa  also  subscanclace  Eqa.  (3)  and  (4)?  Analyses  of  Che  classic  data  re¬ 
ported  In  (6)  and  [7]  (Comp-B  and  Tetryl),  and  [8]  (TNT),  also  substantiated 
Eqs.  (3)  and  (4)  to  a  remarkable  degree.  Tne  more  recent  test  results  for 
PBX-9404  [1]  also  corroborated  Eqs.  (3)  and  (4)  over  a  rather  wide  range  of 
projectile  diameters.  The  constants  (Ci,  C2,  and  C3)  are  different  for  the 
different  explosives  as  shown  In  Tables  1  through  5,  and  Figures  3,  4,  and  5. 

Tables  1-5  list  results  for  Comp-B3,  Comp-B,  ITT,  Tetryl,  and 
PBX-9404,  respectively.  Figures  3  and  4  illustrate  the  variation  of  Pg 
as  a  function  of  Pg.  Figure  5  depicts  Up^^^  tc^/ yA^r  plotted  versus  Pg.  Note 
that  Pg  is  a  function  of  the  Impact  velocity  (Vx)  (Appendix  A).  Pg  was  em¬ 
ployed  as  a  unifying  parameter  because  all  of  the  experimental  data  were  not 
acquired  with  a  single  projectile  material  (Appendix  B).  The  Vx  required  to 
produce  a  certain  Pg  will  be  somewhat  different  for  the  different  projectile 
materials  (mild  steel  and  copper,  or  brass). 

The  values  of  Pg Y^r  Up  t^r/ V^cr  ®  data  were,  of 

course,  sensitive  to  the  basic  Input  §ata  for  Che  explosive  and  projectile 
(Appendices  B  and  C).  Some  examples  which  focused  attention  on  this  fact  are 
as  follows; 

The  initial  computations  for  the  pressed  TNT  data  employed  Cq  and  S 
values  for  cast  TNT  from  [12].  When  the  appropriate  Cq  and  S  values  for 
pressed  TNT  [26]  were  utilized,  the  results  were  much  more  consistent.  Impact 
velocities,  Vx,  for  the  TNT  data  were  difficult  to  determine  precisely  from 
Figure  19  of  [8].  However,  by  making  iterative  computations  with  minute 
changes  in  Vx,  (for  a  given  projectile  diameter),  very  uniform  Pg V^cr  results 
were  obtained.  The  minqte  changes  in  Vx  were  in  the  third  and  fourth  decimal 
places  for  Vx  expressed  in  cm/  -sec.  This  final  set  of  data  for  TNT  (Table  3) 
also  provided  very  uniform  results  for  the  analysis  of  Section  II. B  (Table  8), 
which  is  independent  of 

The  density  of  the  [6]  Tetryl  was  do  ■  1.54  grams/cm^.  However,  the 
density  of  Tetryl  employed  to  acquire  Cx  was  1-655  gram/cm^  [16]  •  This  result 
for  Cx  was  employed  since  Cx  for  Oo  ■  1.54  grams/cm^  was  not  available.  Some 
of  the  Inconsistency  In  the  Tetryl  results  may  be  caused  by  mismatched  Cx 
values.  The  results  for  PBX-9404  appeared  somewhat  more  consistent  when 

Cx  •  50.244  J/cm^  [4]  was  employed  Instead  of  Cx  ■  64.4  J/cm^  from  a  later 
source  [ 14 ] . 

Although  not  thoroughly  Investigated,  the  projectile  shock  and  par¬ 
ticle  velocity  information  (Appendix  B)  Is  obviously  an  Important  Input  also. 
The  steel  projectiles  are  assumed  to  be  low  carbon  mild  steels.  As  such,  they 
may  undergo  a  phase  change  when  shocked  above  about  130  Kllobars  (Kbars).  The 
shock  and  particle  velocity  relationship  will  be  different  above  and  below 
this  transition  if  it  existed  for  the  steel  of  projectiles.  Recommendations 
for  further  Investigation  of  this  possible  steel  projectile  phase  transition 
effect  are  given  In  Section  IV. 


4 


Pros  Che  above  dlecucsion,  and  chat  In  Appendix  C,  it  should  be  obvious 
Chet  obtaining  a  conaistant  set  of  input  information  for  both  the  explosives 
and  the  projectiles  vas  a  primary  concern  of  the  authors.  Some  additional 
methodology  concerns  are  also  expressed  in  the  form  of  action  recoauaenda Cions 
in  Section  IV. 

Additional  relations  involving  Acr  and  the  shock  variables  can  be 
derived  from  Eqs.  (1)  and  (2).  For  instance,  one  such  relation  is: 

Ps_^r  _  ^r(pQ,,  _  Force 

°Pex  Length 

C2^ 

»  —  ■  C4  •  Constant  ,, 

Cl  (5 


Where 

Poex  ~  Initial  unshocked  density  of  the  explosive 
•  Shock  front  velocity  in  Che  explosive 

Tables  1  through  5  list  Pg  Acr/^p  ^  icr  test  data  considered 

herein.  For  certain  explosives,  Che  megnlfude  of  this  parameter  is  very  uni¬ 
form.  The  greatest  deviation  from  uniformity  occurs  with  FBX-9404  for  projec¬ 
tile  C  (19.050).  A  postulated  explanation  of  the  behavior  of  PBX-9404  for  the 
Larger  diameter  projectiles  is  given  in  Section  II. B. 

It  is  worth  noting  that  Eq.  (3)  states  that  Pg  is  inversely  propor¬ 
tional  to V^r *  That  is: 

C2  C2  I 

^  xF  ^ 

This  relation  is  analogous  to  a  similar  expression  for  the  variation 
of  elastic  wave  stress  amplitudes  in  a  slender  conical  rod  [38]  or  to  spheri¬ 
cally  symaatrical  elastic  wave  amplitudes.  In  Section  II. C,  both  conical 
(Eq.  (23))  and  spherical  (Eq.  (25))  volume  relations  are  derived  from  the 
energy  integral  (Eq.  (20)). 
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B.  Projectllt  Crosi-Section  Area  (A)  Rolationa 

In  Appendix  F,  It  was  conjectured  via  analogous  reasoning  that  there 
was  a  relation  between  the  critical  Ignition  Zone  Length  (IZL)  and  the  rtrojec- 
tile  cross-section  dimensions.  This  conjecture  wee  verified  by  empiri..al  num¬ 
erical  experimentation.  The  resulting  expression  (Eq.  (6F)  in  Appendix  F)  is: 

IZL  (U.„  -  Up„)  ter 

- W - 

••  Constant  (6) 

C5  is  a  constant  or  practically  a  constant  for  a  given  explocive. 
ter>  A  have  all  been  previously  defined  in  Section  II .A. 

Table  6  contains  the  computed  values  of  IZL/V^  for  each  explosive 
considered  herein.  These  results  are  graphically  depicted  in  Figures  6  and  7. 
The  average  value  of  IZL/\^  for  each  explosive  is  considered  to  be  C5  for  that 
explosive. 

Equation  (6)  states  that  the  IZL  is  directly  proportional  to  the 
projectile  diameter,  or  more  generally  to  the  square  root  of  the  projectile 
cross-sectional  area  CVX).  Note  that  Eq.  (6)  is  not  based  on  the  critical 
area  (A^r)  concept. 

Regardless  of  its  present  uncertain  physical  interpretation,  Eq.  (6) 
must  define  some  natural  parameter  since  its  excellent  corroboration  for  the 
five  diverse  explosives  considered  herein  can  hardly  be  dismissed  as  fortui¬ 
tous  or  coincidental.  One  of  the  better  correlations  came  from  the  FBX-9404 
data  over  an  extremely  broad  range  of  conditiona.  In  particular,  the  correla¬ 
tion  was  good  for  the  larger  projectiles,  which  exhibited  considerable  devia¬ 
tion  from  the  norm  for  the  critical  area  parameters  (Section  II. A)  and  energy 
input  parameters  (Section  II. C). 

It  is  believed  that  for  the  larger  diameter  projectiles,  PBX-9404 
behaves  in  the  same  manner  as  the  propellant  LOVA-XIA  reported  in  [41].  It 
was  shown  in  [41]  that  ignition  (bum)  and  detonation  Impact  velocities  are 
different  at  the  larger  diameters,  but  are  the  same  at  the  smaller  diameters. 
Also,  the  ignition  onset  always  seemed  to  follow  a  constant  critical  energy 
criteria,  but  detonation  initiation  did  not  do  so  for  the  larger  diameter  pro¬ 
jectiles  (Figure  3,  [41].  See  also  the  discussion  in  [22],  particularly  the 
comments  by  F.  E.  Walker  with  respect  to  possible  low  shock  pressure  effects. 

Other  empirical  relationships  can  be  derived  by  combining  Eq.  6  %rlth 
Eqs.  (1)  through  (5)  of  Section  II. A.  Two  of  the  more  plausible  results  are: 


6 


1 


P,VAcr  *"  A  .  PsVAcr  *  A 

(U*ax  -  ^ct 

-  C2/C5 

-  C6 

■  ConfltanC 

'W  ^Pex  “  W  V  P“Po/  ** 

-  C5/C3 

-  C7 

■  Constant 


(7) 


(8) 


Raaulta  fro«  Eqa.  (7)  and  (8)  ara  tabulatad  In  TabU  6  for  aach  data 
po‘nt.  Thay  art  valid  on  tha  boundary  bttwtan  datonation  and  no-datonatlon. 
Note  th#  raaulta  for  TNT  (Table  6)  where  it  la  believed  tb^t  all  points  are 
on,  or  very  cloae,  to  this  boundary. 

It  w^a  shown  via  Figures  G-1  through  G-5  In  Appendix  G  that 

(P,  -  PjiIq)  dp  -  Constant  (9) 

Where  P^in  d*  the  extrapolated  point  where  1/dp  equals  tero,  or  one- 
dioeaslcnal  conditions  exist.  Since  dp  is  the  projectile  dlaaeter,  Eq.  (9) 
Implies  that: 

(P,  -  P»ln)VA  -  Constant 

-  Cg  CIO) 

If  Eq.  (10)  Is  oultlplled  by  Eq.  (6),  this  yields: 

(Ps  -  ?»ln)  CUs,x  -  0p,x^  -  C5  <  Cg 

(Ps  -  Pain)  CIZL)  -  C9 

-  Constant  (11) 


Nuaerlcal  results,  given  In  Table  7,  substantiated  Eq.  (10)  for  tha 
five  sets  of  data  under  Investigation  (Figures  8  and  9).  Equation  (11)  was 
valid  for  pressed  TNT,  Tetryl  and  PBX-9604  (Figure  10).  CoBp-B3  and  Coap-B 
results  exhibited  a  linear  variation  of  (P»  -  Pain)  (IZL)  as  a  function  of  P#. 
Thus,  for  Coap-B3  and  Coap-B,  C9  of  Eq.  (8)  is  not  a  constant,  but  could  be 
written  as  Cj^q  Cn  *  P*,  wher**  and  Cj^i^  are  constants. 
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Note  Chat  Eqs.  (10)  and  (II)  are  aomawhac  analogoua  to  Eqs.  (3)  and 
(1),  raapactlvaly t  in  Section  II. A.  In  the  saoe  manner,  there  la  a  similar 
analogy  between  Eq.  (6)  in  this  section  and  Eq.  (4)  in  Section  II. A.  These 
equations  are  collected  and  displayed  side  by  side  in  Figure  11. 

C.  Power  and  Energy  Relations 

The  eopirical  retulca  of  Section  II .A,  clearly  demonstrate  the  per¬ 
vasive  Influence  of  the  critical  area  over  a  vide  range  of  conditions  for  a 
given  explosive.  One  basic  relation  (Bq.  (3))  and  two  auxiliary  relations 
(Eqs.  (4)  and  (5))  were  obtained  which  were  functions  of  the  critical  area, 

Che  critical  time,  and  the  shock  variables.  These  relations  were  either  di- 
■ensionless  (Eq.  (4))  or  had  the  units  of  energy  fluence  (Eq.  (1)).  A  basic 
and  persistent  question  was:  How  do  these  relationships  fit  into  overall  gen¬ 
eral  relations  ^ich  govern  the  physical  phenoaenla?  In  other  words,  what  arc 
these  general  relationships? 

Consequently,  in  search  of  increased  physical  insight,  additional 
analysis  and  empirical  numerical  exploration  with  respect  to  the  shock  energy 
and  energy  input  rate  was  performed  as  follows: 


Consider  the  projectile  as  a  piston  with  a  transient  cross-section 
area,  A(c).  A(c)  is  Che  variable  shock  loaded  area  of  the  pro jectlle/explo- 
siva  target  interface  as  described  by  Moulard  ([2],  [3],  and  Appendix  D).  The 
total  area  dependent  high  shock  energy  input  rate  or  power  input  of  the  piston 
at  this  interface  Is: 

f  •  A(C)  P.  Up., 


-  E  -  J/u-sec  (12) 

This  Is  the  abbreviated  form,  coamensurate  with  the  critical  energy 
fluence  criteria  as  commonly  employed  (Eq.  C-1  of  Appendix  C).  The  exact 
expression  consMnsuraCe  with  Eq.  C-3  in  Appendix  C,  is: 


dE 

It 


A(t) 


-  A  (t) - : - 

®Oez  (^’■ex  -  ^Pex) 


(13) 
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In  Che  following  analysis,  Eq.  (12)  will  be  employed  since  ic  is  com¬ 
patible  with  critical  energy  flux  density  (Pg  tcj.  -  C^)  which  is  utilised 
extensively' 


chat : 


For  simplicity,  consider  a  cylindrical  projectile  of  Radius,  R,  such 
A(C)  -  w  r(t)2 


-  IT  (r2  -  2R  Ugp  t  +  (Ugp  t)2)  (14) 


Since 

r'vt)  -  R  -  Ugp  t 

(15) 

At  an  arbitrary  time,  t]^,  the  energy  rate  is: 

f  ■  A(ti)  P.  -  (i),j 

(16) 

(E)(^  signifies  the  tisM  derivative  of  E,  evaluated  at  C^.  When 
c^  ■  0,  then; 

(E)to  -  A  Ps  Vx 

(E)(^  is  Che  oaxlBum  input  energy  rate.  'Its  nuaerical  values  are 
listed  in  Table  8  for  the  five  sets  of  data  under  investigation. 

It  Is  noteworthy  that  whan  •^1  •  tcr»  chan  (E)tgf  can  be  axpraasad  as 
follows  via  the  empirical  relations  astabllshad  lu  Section  II. A: 


dE 

St 


“cr 


**•  ^Pex 


Acr 
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gg^j 


Const . 


■  (E)tcr 

Thus,  froB  Eq.  (18],  (E)tgj.  Inversely  proportional  to  the  explo¬ 
sive  shock,  velocity;  so  it  is  not  a  conatant.  Howtvar,  it  should  not  exhibit 
order  of  nagnitude  variations  as  delineated  belov.  The  iaportanca  of  the 
shock  velocity  influence  is  indicated  by  the  following  postulated  upper  and 

lower  bounds  for  (i)f 

'•cr 


(E)tcr  - 


Where 

>  Blniaua  shock  velocity  required  for  self-sustained  phonon 
fission  or  crystal  lattice  breakup  and  disintegration  (Appendix 
G,  (46]  and  (47)). 

■  explosive  detonation  shock  velocity  (Appendix  F,  Table 

F-1). 

These  sstiBStes  for  the  limits  on  dependent  on  the  assump¬ 

tion  that  C2  remains  constant  at  these  extreme  limits. 

e 

Computed  values  for  the  shock  energy  derivative  pre¬ 

dicted  upper  sod  lower  bounds  (Eq.  (19)  and  Table  G-2)  are  listed  in  Table  8 
and  shown  versus  Pg  in  Figures  12,  13,  and  14  for  all  five  sets  of  experimen¬ 
tal  data.  bounded  reasonably  well  by  Eq.  (19)  with  the  exception  of 

PBX-9404  at  its  upper  limit.  See  the  comments  in  Section  II. B  with  respect  to 
PBX-9404  and  the  large  diameter  projectiles. 

f 

The  faired  values  for  Comp-B3  (Figure  12)  and  Tetryl  (Figure 

14)  are  somewhat  speculative  because  of  the  scatter  in  the  computed  results. 
More  refined  computations  are  planned  as  described  in  Section  IV. 

s 

After  observing  (E)tj.j.  versus  P,  in  Figures  12,  13,  and  14,  Dr.  Joe 
Foster  [42]  remarked  that  the  trend  was  similar  to  published  results  for 
gaseous  explosive  detonation  power  versus  energy  criteria,  since  Pg  Is  energy 
per  unit  volume.  Susbequent  inquiries  ravealed  (43],  (44],  and  (45]  which 
delineate  power  and  energy  requirements  for  gaseous  explosive  detonation.  The 
similarity  of  Figure  6  in  [43]  for  gaseous  detonation,  and  Figures  12,  13,  and 
14  in  this  report  for  solid  explosive  detonation,  illustrates  rather  vividly 
the  analogous  energy  and  power  requirements  for  both  gaseous  and  solid  explo¬ 
sive  detonation. 


C2^ 


Ooex  “"f 


as  t(.f 


Poex  ^“ex 


*■  ^cr  +  0 


(19) 
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That  Is,  for  either  solid  or  gaseous  explosives  subjected'to  low 
shock  (or  energy)  situations,  a  high  energy  rate  or  power  input  is  required  to 
initiate  irreversible  reactions  leading  to  detonation.  Conversely,  for  both 
solid  and  gaseous  explosives,  high  shock,  pressure  (or  energy  input)  require  a 
relatively  low  power  input.  In  any  case,  for  a  given  shock  pressure  loading 
(or  energy  input)  at  least  a  certain  lainiaum  energy  input  rate  (or  power)  laust 
be  applied  to  either  gaseous  or  solid  explosives. 

From  Eqs.  (12)  and  (lA)  the  total  high  shock  energy  (J)  may  be  cob- 
puted  for  a  cylindrical  projectile  of  radius,  R.  This  is  accomplished  as 
follows  to  find  E  (ti)  •  El  for  an  arbitrary  time,  ti: 

r^l 

El  .  J  A(t)  P.  Up„  « 

0 

/I 

-  J  (r2  -  2  R  U,p  t  +  U.p2  t2)  dt 

o 


-  (r2  +  R  Ri  +  Ri2)  (Up^^  ti)  P, 


Where  Ri  -  R  -  U,  tj.  (21) 

^  • 

In  geaeral,  amy  b«  vriCtcn  as: 

El  -  ^  (A  +Va  *  Ai  Ai)  (Up^^  ti)  P.  (22) 

Where 

A  ■  projectile  cross-section  area 

Ai  ■  high  shock  loaded  interface  area  at  tlise  ti  (Appendix  D) 
Equations  20  and  22  have  an  interesting  interpretation.  The 

quantity: 

Voli  -  y  (A  VA.  Ai  +  Ai)  (Up^^  ti) 

"  ^\l^  ^^Pex  <23) 

is  the  volume  of  a  general  conical  frustum  whose  length  is  Up^^  ti 
and  whose  end  areas  are  A  and  Ai  (Figure  15).  Ay^  is  the  volume  weighted  area 
of  the  frustrvm  cross-section.  For  a  cylindrical  projectile: 

ti  -  (R  -  Ri)/U8p 

where  Ug^  •  shock  velocity  in  the  projectile  (24) 
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So  that : 


Voli 


i  (r2  +  R  *  Rl  +  R]_2)  (r  -  R^) 


'Pex 


-  (r2  -  Ri3) 


^Pex 


(25) 


In  this  form,  Volj^  la  ona-fourth  of  the  difference  between  two 
epherical  volumes  cIms  the  velocity  ratio. 

The  average  energy  Input  rate  from  t  -  o  to  t  ■  t^  Is: 


lEi) 


av 


ti 


■  Av]^ 

Pa 

"Pex 

When 

ti  -  0, 

Eo 

■  0  and 

dE 

-  A  Pg 

0 

^Pex 

If  ti 

■  tcr» 

then: 

Ecr  • 

(A  + 

■VT 

•  Acr 

A^r)  (ypgj  ter)  Pg 

■ 

(Avcr) 

Es 

"Pex 

ter 

- 

(Av,,) 

Cl 

-  (Volcr)  Ps 

-  1  (r3  -  R^.j.3) 

3  U.p 

-  Y  -  Rcr^)  POp  Upp  Up„ 


(26) 


(27) 


(28) 


The  average  critical  energy  Input  rate  Is: 


(Ecr) 


av 


'cr 

-cr 


■cr 


-  (Av,,)  P. 


'Pax 


(29) 
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If  tj  -  T  *■  R/L'gp,  t(T)  ■  0  and  Aq.  -  0.  Thus: 
tr  -  i  Uh  (Up„  T)  t, 

■  7  <“P.. 


■  (“Pex 

-  (Vol^)  P, 

•  (■  ■•) 


(POp  l^Pp  Up.x) 


OO) 


This  ralatlOQ  has  both  conical  and  spherical  volune  connotations. 
The  average  Input  rata  for  E/p  Is: 

«  E<t> 

<^)av  -  T 


T 


'Pex 


(31) 


Table  9  lists  computation  results  for  (Eq.  (28))  for  the  five 

sets  of  data  being  examined.  was  employed  to  compute  Ecr  (▼!»  Eq.  (28)) 

and  (Ecr)  (via  Eq.  (29))  and  these  quantities  are  given  In  Table  9.  Table  9 
also  lists  numerical  results  for  E.J.  (via  Eq.  (30))  and  (E^)^^  (via  Eq.  (31)). 


What  Is  the  "appropriate  energy  input  rate  7"  Is  It  Acr  Ps  Up^jj  or 
Ayer  Pg  Up^^?  The  average  value,  Ay^r  ^s  ^P*x’  Interval,  tert 

is  larger  than  A^r  P»  ^p.x’  *  more  conservetlve  measure  of  the  power 

input  requirements.  Certainly  A^r  Ps  Up.x  ^  considered  an  absolute 

bare  nlnlDum  requirement. 

Consider  the  Impact  geometry  depicted  In  Figure  16  for  typical  tests 
where  the  areas,  A^Of  of  both  the  flat  plate  projectile  and  explosive  target 
are  ouch  greater  than  A^r  ^or  the  particular  test  conditions.  Figure  16. A 
depicts  conditions  at  the  Instant  of  contact  (t-o). 
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In  Figure  16. B,  the  time  after  Impact  la  t^j  ■  ZW/Ug^  ao  the  ahock  In 
the  projectile  has  reflected  back  and  Is  at  the  projectile  target  Interface. 
Also,  at  the  edges  of  the  projectile/target  interface,  rarefaction  waves  have 
formed,  similar  to  the  small  projectile  Impact  tests  (Figure  1).  The  high 
pressure  area  at  the  Interface  Is  tRi^,  or  Ai  which  Is  much  larger  than  A^^  or 
even  for  that  particular  pressure  level.  In  fact  A]^  Is  approximately 

equal  to  Aqq  and  the  critical  area  Is  greatly  exceeded.  Thus,  the  Ignition 
phenoMoa  is  essentially  area  independent.  That  is,  if  Eq.  (22)  is  applied  to 
this  condition,  then: 

Ecr  -  y  (A,  +Va.  *  Ai  +  Ai>  P,  ter 

*  ^  “P,x  ^cr  (32) 


Additional  tests  with  large  specimen  and  different  place  thickness 
(t(.r  ■  ZV/Ug  )  yield  the  familiar  constant  energy  per  unit  area  criteria  (for 
certain  explosives,  or  course). 


Es  (Jpex 


••  Constant 


-  Cl 


(33) 


This  is  area  indepandenc  and  is  a  special  case  of  the  more  general 
area  dependent  pover/energy  concept  postulated  in  this  section. 
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III.  CONCLUSIONS 


This  investigacioa  has  revealed  simple  empirical  relaclonshlpa  becween 
^he  crosa-sectlon  section  area  (A)  and/or  the  critical  cross-section  area 
(A^r)  of  a  siaall  flat  faced  projectile  and  certain  quantities  (Pg,  Up^^,  t^j.) 
which  are  critical  Impact  shock  variables  for  fWe  different  explosives.  The 
uniformity  and  consistency  of  the  prltaary  empirical  relations  (Eqs.  (3),  (4), 
(6),  (10),  and  (11)}  has  generated  considerable  confidence  In  their  validity, 
although  their  general  physical  meaning  may  not  be  clear.  It  Is  believed  that 
they  are  portions  of  overall  general  relationships,  which  are  physically 
Bieanlngful.  For  Instance,  Eq.  (3)  (Pg  ■  C2)  was  employed  to  simplify 

e 

Eq.  (18)  for  (2)tj,j,  *0  that  the  Influence  of  the  shock  velocity  was  clearly 
delineated. 

The  analysis  of  the  area  dependent  shock  energy  Input  rate  (or  power) 
strongly  suggests  that  It  Is  a  fundamental  parameter.  For  a  given  shock 

e 

pressure  Input,  at  least  a  minimum  energy  rate  must  be  applied  for  a 

length  of  time  equal  to  t^^.  This  detonation  criterion  for  solid  explosives 
Is  analogous  to  results  reported  In  the  literature  [43],  [44]  and  [4S]  con¬ 
cerning  gaseous  explosive  detonation.  It  Is  shown  also  that  the  area  depen¬ 
dent  power  Input  criteria,  and  the  constant  energy  (Pg  t^r)  criteria  are 

Indistinguishable  when  true  one-dlaenslonal  conditions  exist  (A  »  Acr)  • 


One  Important  spin-off  which  evolved  from  this  Investigation  has  been 
the  application  of  a  crystal  lattice  fracture/dlalntegratlon  criteria  [46]  to 
the  prediction  of  minimum  shock  pressure  required  to  detonate  an  explosive. 
Some  pertinent  results  are  contained  herein  (Section  II. C  and  Appendix  G),  but 
the  analysis  and  concept  are  being  published  In  a  separate  report  [47].  Cer¬ 
tain  phenomena  observed  In  shocked  Inert  materials  may  also  be  explained  via 
the  crystal  lattice  disintegration  criteria. 
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IV.  RECOMMENDATIONS 


A  coQglderable  ai&ouQC  of  conputacional  ref Ineiaeat':  and  numerical  axperi- 
mantation  remain  Co  be  accompliahed. 

a.  Firac  of  all,  Che  influence  of  a  possible  polymorphic  phase  cransfor- 
macion  efface  in  che  mild  steel  projecciles  ac  approximaCely  130  Kbars  shock 
preasure  needs  Co  be  investigaced.  This  rn^y  ..‘aquire  chat  cvo  different  func- 
Cions  for  the  steel  shock-particle  velocity  relationship  be  employed  (one 
above  130  Rbers  and  one  below  130  Rbers).  There  are  a  number  of  papers  and 
reports  on  this  effect  (not  cited  herein)  which  must  be  consulted*  The  preo- 
ent  analysis  employed  a  shock-particle  velocity  function  which  is  best  at 
somewhat  higher  shock  pressures  than  shown  for  any  of  the  data  considered 
herein  (above  200  Kbars).  Some  of  che  nonunif ormicy  in  the  present  results 
for  steel  impacts  may  be  caused  by  an  incorrect  shock-particle  velocity  rela¬ 
tion  rather  than  any  critical  energy  criteria  limitations  (particularly  for 
the  lower  shock  pressures  and  larger  projecciles).  The  criterion  for  this 
phase  change  analysis  will  be  a  comparison  of  the  uniformity  and  consistence 
of  results  with  those  presently  obtained. 

b.  Determine  It  possible,  via  literature  search,  etc.,  the  appropriate 
values  or  relationships  for  the  rarefaction  shock  velocity  for  both  the  ex¬ 
plosive  target  and  projectile  materials.  For  example,  [37]  contains  such  a 
relationship  in  terms  of  the  shock  and  particle  velocities.  This,  and  ocher 
feasible  relations  or  values,  should  be  employed  in  the  following  computa¬ 
tions: 

1.  Recompute  A^^  and  associated  relations  using  Che  rarefaction 
shock  velocity  for  che  projectile  material. 

2.  Recompute  A^^^  and  associated  relations  using  che  rarefaction 
shock  velocity  for  the  explosive  target  material. 

3.  Recompute  A^^  and  associated  relations  using  the  compression 
shock  velocity  of  che  explosive  target  material. 

Compare  che  paragraphs  I,  2,  and  3  results  with  Che  present  results 
(which  used  the  projectile  compression  shock  velocity)  for  uniformity  and  con¬ 
sistency.  The  object  of  all  this  numerical  and  input  experimentation  is  to 
conform  more  closely  Co  actual  physical  conditions.  This  should  yield  more 
uniform  results  from  Eqs.  (3),  (4),  (6),  (10),  and  (11). 

Note  that  computation  of  via  che  rarefaction  shock  velocity  in  the 
projectile  is  suggested  by  Moulard  [2].  However,  employing  the  shock  velocity 
of  the  explosive  material  (2  or  3,  above)  may  be  physically  more  realistic 
since  che  explosive  also  unloads  at  Che  projectile  interface  via  a  rarefaction 
wave  (Figure  1).  Thus,  che  explosive  rarefaction  velocity  may  be  SK>re 
appropriate  for  computing  A<.f  than  the  projectile  rarefaction  shock  velocity. 

The  application  of  Fitzgerald's  crystal  lattice  fracture  criteria  [46]. 
to  shocked  induced  explosive  detonation,  was  deemed  important  enough  Co  war¬ 
rant  documentation  in  a  separate  report  [47],  which  contains  several  sugges¬ 
tions  for  future  work. 
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Figure  2.  Impact  velocity  required  co  initiate  detonation  In  five  bare  explosives. 
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Figure  3.  r-g  \  f^cr  variation  for  CoBp-B3,  Coap-B  and  PBX-9404  . 
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Figure  6.  IZL/yA  variation  for  Cottp-B!i ,  Comp-B^  and  TNT. 
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figure  7.  IZV/yT^  variation  for  Tetryl  and  PBX-9404 ■ 
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Figure  10.  (Pa  -  Pmln)  (IZL)  variattoo  for  Co«p-B3,  Comp-B,  TNT,  Tetryl,  and  PBX-940A . 
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Figure  11.  Analogous  relations. 
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itatlons.  (Sketch  not  to  scale.) 
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Figure  16.  Impact  conditions  for  classic  Pa  t^,^  tests. 
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TABLE  1.  Critical  Area  Coaputatlona  for  Coaip-B3  .  (Concluded) 
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TABLE  5.  Critical  Area  Coaputatlona  for  PBX-9404. 
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TABLE  5.  Critical  Area  Computations  for  PBX-9404  .  (Concluded) 
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Note:  The  averages  for  C2,  C.j,  and  do  not  Include 

results  for  projectile  C-1.9050. 
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TABLf  9.  Tabulation  of  Ayer  Coaputatlona . 
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TABLE  10.  Energy  and  Average  Energy  Rates  for  t^.,.  and 


^  \ 

«l  • 

O  r? 


M 

®  s « 

•-  12  * 
l<  \n  ^ 


• 

j  u  ^ 

d  u  9 

I  P  O 


li- 


SOOOOO  OOC  PQOQ 

o5000  OOO  0©©0  <N4«Nfsr^rg<N 

OOOOOO  <N<M«S<MfS  fNfSCNCN  QQPPPP 

CDOOflO 


^  ^ 

O  (N  r«<»  0 


40^<s  fM««9vtAkA 


rs<«^c>4 


8i0^fN.r^<n  ^  ^  m 


r^lA® 


000^0^  OOOO^ 


^  fn  sc  ^  ^  ^  fs®*^ 

r4  r^  ^  ^  ^  f^^®  ®®<«fs^ 


O  ©  ^  tA 
®  5  r>  ^ 


9^Mtnool^T 

®®^W^  «C0«"®^0 

r>9^®^  0m  ^  ^  ^  C<4 


y\  ir\  m*  ^  ^ 

O  9*  ^  1^  ^ 


2<^®^^9i  0^  m  ^^®irk»A 

Or^O^soi  O®^ 


»r\  ^  r^ 

^  ® 


®l*^®  ®®rs4r«k  rvh%®p«»*n® 

®^0  ^w^®®  tr\0*^®0® 


O  ^ 

'^2^2 


00  ^  r^  rm  ^O®® 

•M  ^4S®®0  «-«fSm 


o  ^  ®  ® 


(M®f^«s<s®  ©r»^  ^  ^  m  ^  o  O  ^  00  ^  ^  ^  ^ 

r^®OiNr<«?»  iN^f^  ^rsf^O 

C^viMf-MfNfN  — tn  ^  ^  tr  ^  ir>  ^  OOO^c^ 

OOOOOO  OOO  OOO*^*^  OOO^  OOOOO 


®®©—  5^®0  tn  ^  r4  ^  IT  is®i<r®  in— i^-*® 


fN®9^9t^9t  «^®®  ©Ck  — 

<n— ^  mm  ^ 
^  tn  tr>  rm,  rm  csr^^ 


®  f>*>  —  9^ 
in  m  O  04 


^  /"s  ^— »©©  — fc— .p 

— «  — »0^  in©©^  ©®®0® 

p«^  ^  OO©  r^o#r^9i  ^f<fcino4  — 

—  —  ^  ©w^***  pK®**  f*4^in®» 

n%MMOM^  *04^  ••©0410  ••04®  ••••© 

®  —  <nr*  —  —  r>r**  — 

W'  w  w  W  W  W  W  V  W  W  W'  W  W  ^  W  ^  W  w 

0«<U^0  OOU  OOUOU  OOOU  OUUUO 


T  T  T  7  I  I  1(1 

%^rSx%  x^z  ttttt 

883888  888 


3S2333 

gggg  assssa 


L 


ZII  55$*-0  OC-KZ  (050*61)3 


REFEtlENCES 


1.  3ahl,  K.  L*,  Vantloe,  H.  C.,  and  Walngart,  R.  C.,  "The  Shock  InlClatlon 
of  Bare  and  Covered  Explosives  by  Projectile  Impact,"  paper  in  the  pro** 
ccedings  of  the  Seventh  Symposium  (International)  on  Detonation, 

16-19  June  1981,  proceedings  published  as  NSWC/HP-82-334,  pp.  325-335. 

2.  Moulard,  Henry,  "Critical  Area  Concept  for  Initiation  of  a  Solid  High 
Explosive  by  the  Impact  of  Small  Projectiles."  Paper  presented  at  the 
Seventh  International  Colloquium  on  Gas  Dynamics  of  Explosions  and 
Reactive  Systems,  PRG,  August  1979,  published  in  the  AIAA  Progress  in 
Astronautics  and  Aeronautics.  Vol.  75,  Gas  Dynamics  of  Detonations  and 
Explosions,  1981,  pp.  196-302. 

3.  Koulard,  Henry,  "Critical  Conditions  for  Shock  Initiation  of  Detonation 
by  Small  Projectile  Impact,"  paper  in  the  proceedings  of  the  Seventh 
Symposium  (International)  on  Detonation,  16-19  June  1931,  proceedings 
published  as  NSWC/MP-82-334,  pp.  316-324. 

4.  Walker,  P.  E.  and  Wesley,  R.  J.,  "Critical  Energy  for  Shock  Initiation  of 
Heterogeneous  Explosives,"  Explosivstoffe,  Nr.  1/1969,  pp.  9-13. 

5.  Glttings,  E.  F.,  "Initiation  of  a  Solid  Explosive  by  a  Short  Duration 
Shock,"  paper  in  the  proceedings  of  the  Fourth  Symposium  (International) 
on  Detonation,  12-15  October  1965,  published  as  report  ACR-126,  pp.  373- 
380. 

6.  Slade,  D.  C.  and  Devey,  J.,  "High  Order  Initiation  of  Two  Military 
Explosives  by  Projectile  Impact,”  U.S.  Army  Ballistic  Research  Labora¬ 
tory,  Aberdeen  Proving  Ground,  HD,  Report  No.  1021,  July  1957. 

7.  Weiss,  Milton  L.  aud  Litchfield,  Elton  L. ,  "Projectile  Impact  Initiation 
of  Condensed  Explosives,”  U.  S.  Department  of  the  Interior,  Bureau  of 
Mines,  Report  of  Investigations  6986,  1966.  (Contains  [6]  data.) 

8.  Eldh,  D.,  Person,  B.,  Ohlin,  B.,  Johansson,  C.  U.,  Ljungberg,  S.,  and 
SJolin,  T.,  "Shooting  Test  with  Plane  Impact  Surface  for  Determining  the 
Sensitivity  of  Explosives,”  Explosivstoffe ,  Nr.  5/1963,  pp.  97-103. 

9.  Roth,  J.  P.,  "The  Initiation  of  Explosives  Caused  by  Action  Mechanism,'’ 
paper  presented  at  the  Third  Symposium  on  Chemical  Problems  Connected 
With  the  Stability  of  Explosives,  Ystad,  Sweden,  28-30  May  1973. 

10.  Walsh,  J.  H.,  Rice,  M.  H. ,  McQueen,  R.  G. ,  and  Yarger,  P.  L.  ,  "Shock 
Wave  Compression  of  Twenty-Seven  Metals.  Equation  of  State  of  Metals," 
Physical  Review,  Vol.  108,  No.  2,  15  October  1957. 

11.  McQueen,  R.  G.  and  Marsh,  S.  P.,  "Equation  of  State  of  Nineteen  Metallic 
Elements  from  Shock  Wave  Maesurements  to  Two  Megabars,"  Journal  of  Applied 
Physics,  Vol.  31,  No.  7,  July  1960. 


49 


REF^'.RENCES  (Continued) 


12.  Jaaeson,  k.  L. ,  Boyls,  V.  M.,  aud  Sulcanoff,  M. ,  ’'DeCenalnation  Of  Shock 
UugonioCa  for  Several  Condensed  These  Explosives,"  presented  at  the  Fourth 
Symposlun  (International)  on  Detonation,  IZ-'IS  October  1963,  proceedings 
published  as  ACR-126,  pp.  241-247. 

13.  Llndstroa,  I.  E.,  "Planar  Shock  Initiation  ot  Porous  Tetryl,"  Journal 
of  Applied  Physics,  Voi.  41,  Nuaber  1,  January  19'^0,  pp.  337-350. 

14.  Groan,  L.  G.,  Nldlck,  Jr.,  R.  J.,  and  Walker,  F.  E.,  "Critical  Shock 
InltlatloD  Energy  of  PBE-9404,  A  Nee  Approach,"  Lawrence  Llvcraore 
Laboratory,  University  of  California,  LIveraore,  CA,  Report  UCRL-51322, 

25  January  1974. 

15.  Dobratz,  3.  H. ,  "LLNL  Explosives  Handbook,  Properties  of  Chealcal  Explo¬ 
sives  and  Explosive  Slmilants,"  Report  UCRL-52997,  16  March  1981. 

Lawranca  Llveraora  Laboratory,  University  of  California,  Llveroore,  CA, 
94550. 

16.  Walker,  F.  B.,  Green,  L.  G.,  and  Nldlck,  Jr,  S.  J.,  "Critical  Energy  for 
Shock  Initiation  of  Tetryl  And  A-5,"  Lawrence  Llveraore  Laboratory 
University  of  California,  Llveraore,  CA,  Report  UCID-16469,  14  March 
1974. 

17.  Roalund,  L.  A.,  Watt,  J.  W.  and  Colebum,  N.  L. ,  "Initiation  of  Warhead 
Explosives  by  tbe.Iapact  of  Controlled  Fragaents  I.  Noraal  lapact," 
NOLTR-73,  12h,  15  August  1973,  Naval  Ordnance  Laboratory,  White  Oak 
Silver  Spring,  MD,  20910. 

18.  Howe,  P.,  Frey,  R. ,  Taylor,  B.,  and  Boyle,  V.,  "Shock  Initiation  and  the 
Critical  Energy  Concept,”  paper  presented  at  the  Sixth  Syapoelua  (Inter¬ 
national)  on  Detonation,  24-27  August  1976,  proceedings  published  as 
ONR-ACR-221,  pp.  11-14. 

19.  Titov,  V.  H. ,  Sllvestrov,  V.  V.,  Kravtsov,  V.  V.,  and  Stadnltshenko,  I*  A., 
"Investigation  of  Soac  Cast  TNT  Properties  at  Low  Teaparatures,"  IBID, 

pp.  36-46. 

20.  Stresau,  R.  E.  and  Kennedy,  J.  S.,  "Critical  Conditions  for  Shock 
Initiation  of  Detonation  In  Real  Systeaa,"  IBID,  pp.  68-75. 

21.  Hayes,  D.  B.,  "A  P°T  Detonation  Criterion  froa  Tharaal  Explosion  Theory," 
IBID,  pp.  760-81. 

22.  "General  Discussion  on  Shock  Initiation  and  P^/T,"  1.  Coaaents  by  F.  E. 
Walker;  2.  CoaoMnts  by  M.  Cowperthwalte;  3.  Coaaents  by  R.  Frey  and 

P.  Howe,  IBID,  pp.  82-88. 

23.  deLonguevllle ,  Y.,  Fauqulguon,  C.,  and  Moulard,  H.,  "Initiation  of 
Several  Condensed  Explosives  by  a  Given  Duration  Shock  Wave,"  IBID, 
pp.  105-114. 


50 


REFERENCES  (Continued) 


24.  Uoaodel,  C.A.,  Hunphre/t  J<S*.  Weingarc,  R.  C. ,  Lee,  R.  S.,  and  Kraoter,  P. 
’’Shock  Initiation  of  TATB  Fonaula clone,"  paper  praeented  at  the  Seventh 
STopoclun  (Intaixatlooel)  on  Decooetioa,  16-19  June  1981,  publiehed  ae 
NSWC-rtP-82-334,  pp.  425-434. 

25.  Khewainov,  B.  A.,  Borisov,  A.  A.,  Bnsolaav,  B.  S.,  and  Kurotkov,  A.  I., 
"Two-Pha»'e  Vtsco-Plaatic  Model  of  Shock  Initiation  of  Detonation  in  High 
Deneity  Precsed  Explosives,'  IBID,  pp.  435-447. 

26.  ilaasey,  J.  B.  and  Popolato,  A.,  'Analysis  of  Shock  Wave  and  Initiation 
Data  for  Solid  Explosives,"  paper  presented  at  the  Fourth  Syapoalua 
(International)  on  Detonation,  12-15  October  1965,  published  as  ACR-126, 
pp.  233-238,  by  the  Office  of  Naval  Research,  Department  of  the  Navy, 
Arlington,  VA. 

27.  iCictlakowaky ,  C.  B. ,  'Initiation  of  Detonation  in  Explosives,'  Third 
Symposium  oc  Combustion  and  Flame  and  Explosion  Phenomena,  published 
by  the  Williams  and  Wilkins  Company,  Baltimore,  MD,  1949,  pp.  560-565. 

28.  Ubbelohne,  A.  R.,  "Transition  from  Deflagration  to  Detonation:  The 
Physico-Chemical  Aspects  of  Stable  Detonation."  IBID,  pp.  566-571. 

29.  Pandow,  Mary  L.,  Ockert,  K.  F.,  and  Shuey,  H.  M.,  'Studies  of  the 
Diameter  Dependence  of  Detonation  Velocity  in  Solid  ConposlCa 
Propellant.;;  I.  Attempts  to  Calculate  Reaction  Zone  Thickness,' 
paper  presented  at  the  Fourth  S3rmposium  (International)  on  Detona¬ 
tion,  12-15  October  1965,  pp.  96-101. 

30.  Pandow,  Nary  L.,  Ockert,  K.  F,  and  Pratt,  T.  U. ,  "Studies  of  the 
Diameter  Dependence  of  Detonation  Velocity  in  Solid  Composite 
Propellants;  II<  Prediction  of  Failure  Diameter,'  paper  presented 
at  Che  Fourth  Symposium  (International)  on  DaConaCion,  12-15  October 
1965,  pp.  102-106. 

31.  Gordon,  W.  E.,  "Detonation  Limits  in  Condensed  Explosives,"  paper 
presented  at  the  Fourth  Symposium  (XotemaCional)  on  Detonation, 

12-15  October  1965,  pp.  179-197. 

32.  Kusakhbe,  Hasao  and  Fujlvara,  Shuxo,  "Effects  of  Liquid  Diluents 
on  Detonation  Propagation  in  Nicrome thane,"  paper  presented  in  Che 
Sixth  Symposium  (International)  on  Detonation,"  24-27  August  1976, 
pp.  133-142,  published  as  Report  ACR-221,  by  the  Office  of  Naval 
Research,  Departmant  of  the  Navy,  Arlington,  VA. 

33.  Campbell,  A.  W.  and  Engelka,  Ray,  "The  Diameter  Effect  in  High- 
Density  Heterogeneous  Explosives,"  paper  presented  in  the  Sixth 
Symposium  ( Intemational)  on  Detonation,  24-27  August  1976, 

pp.  642-652. 


51 


REFERENCES  (Continued) 


34.  Hardy,  J.  R.,  Karo,  A.  H.  and  Ualkar,  F.  E.  (1979),  'The  Molecular 
DToanics  of  Shock  and  Detonation  Phenomena  In  Condenaed  Matter." 
Gaadynamlca  of  Detonatlona  and  Exploelona;  AIAA  Progresa  In  Aatro- 
oautlca  and  Aaronautlca,  (edited  by  Bowen,  Manson,  Oppenhela,  and 
Soloukhln),  Vol.  pp.  209-225,  AIAA,  New  York. 

35.  Rardy,  J.  R. ,  Karo,  A.  M.  Walker,  F.  E.,  and  Cunnlnghaa,  V.  G. ,  (1981), 
"The  Study  Of  Shocked-lnduced  Signals  And  Coherent  Effects  In  Solids 

by  Molecular  Dynamics.”  GaadynarJ.cs  of  Detonations  and  Explosions; 

AIAA  Progress  in  Astronautics  nnd  Aaronautlca,  (edited  by  Bowen, 

Manson,  Oppenhela,  and  Soloukhln),  Vol.  87,  pp.  9-21,  AIAA,  New  York. 

36.  Cole,  Robert  8.,  Underwater  Explosions,  Dover  Publications,  Inc., 

180  Varlck  Street,  New  York,  NY,  1965.  pp.  20-30,  145,  233-234, 

First  published  by  Princeton  University  Press  In  1948. 

37.  Frey,  R. ,  Melanl,  G.,  Chawla,  H. ,  and  Triable,  J>,  "Initiation  of 
Violent  Reaction  by  Projectile  lapact."  paper  presented  at  Che 
Sixth  Syaposlua  (International)  on  Detonation,  24-27  August  1976. 
pp.  325-335.  Published  as  ACR-221  by  the  Office  of  Naval  Research, 
Departaent  of  the  Navy,  Arlington,  VA. 

38.  Laodon,  J.  W.,  and  Quinney,  H.,  "Experlaents  with  the  Hopklnson  Pressure 
Bar,"  Proceedings  of  the  Royal  Society,  A,  Vol.  103,  1923,  pp.  622-643. 

39.  Green,  Leroy,  "Shock  Initiation  of  Explosives  by  the  lapact  of  Saall 
OlasMter  Projectiles,"  paper  In  the  Seventh  Syaposlua  (International) 
on  Detonation,  16-19  June  1981,  pp.  273-277.  Published  as  N8WC-HF- 
32-334,  Naval  Surface  Weapons  Center,  White  Oek,  Silver  Spring,  MD. 

40.  Walker,  F.  M.,  Karo,  A.  M.,  and  Hardy,  J-  R.,  "Coaparlson  of  Molecular 
Dynaaics  Calculations  with  Observed  Initiation  Phenoaeoa."  IBID, 

pp.  777-788. 

41.  Andersen,  W.  E. ,  "Projectile  lapact  Ignition  Characterlstlce  of 
Propellants,"  Report  3^20-F,  prepared  by  Shock  Hydrocynaalcs  Division, 
Whittaker  Corp.,  4716  Vineland  Avenue,  North  Hollywood,  CA,  91602,  for 
Che  U.S.  Aray  Research  Office,  P.  0.  Box  12211,  Research  Triangle  Park, 
NC,  27709,  Hay  1980,  (AD  Nuaber  Is  A083234). 

42.  Vl^lC  and  discussion  with  Dr.  Joe  Foster  (AD/ALJW),  Eglln  Air  Force 
Base,  Florida,  4  February  1988.  Subsequent  telecon  Inquiries  were 
oade  to  Professors  R.  A.  Strehlow  (University  Of  Illinois,  Urbaua,  IL) 
and  J,  H.  Lee  (McGill  University,  Montreal,  Canada). 
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for  Direct  Initiation  of  Unconflned  Gaseous  Detonations,"  paper  In  Che 
Fifteenth  Syaposlua  (lacernational)  on  conbustlon,  Toshl  Center  Hall, 
Tokyo,  Japan,  25-31  August  1974,  pp.  53-67  of  the  papers  organized  and 
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APPENDIX  A 

REUTIONS  BETWEEN  SHOCK  PRESSURE  (P,),  PARTICLE  VELOCITY  (Up), 
SHOCK  VELOCITY  (U,),  AND  PROJECTILE  IMPACT  VELOCITY  (Vj) 


Consider  the  planar  Impact  geometry  depicted  In  Figure  A-1  near  the  cen¬ 
ter  of  the  projectile.  The  shock  pressure  (Pg)  Is  a  function  of  the  material 
density  (Po)f  shock  velocity  (Ug),  and  particle  velocity  (Up). 

Pg  -  Po  -  Po  U,  Up  sP.  (A-1) 

Experimental  data  reveal  that  the  shock  velocity  Is  a  linear  function  of 
the  particle  velocity  for  many  materials.  Cq  and  S  have  been  determined  for 
many  metals,  plastics  and  explosives.  (Appendix  B.)  Cg  Is  generally  the  bulk 
modulus  elastic  wave  velocity. 


U,  -  Co  +  SUp  (A-2) 

The  projectile-target  Interface  velocity  Is: 

Vl  -  Vi  -  Upp  -  Up^  (A-3) 


Interface  contact  requires  that  the 
projectile  be  equal. 

^•t  -  ^-p 

P*^  “  Po^  ^Pt 

^Sp  ■  POp  Upp  (Cp  +  Sp  Upp) 
Equations  (A-3),  (A-4),  (A-5),  and  (A-6) 

-B  +  Vb2  -  4AC 

V - nj — 


shock  pressure  In  both  target  and 

(A-4) 

(A-5) 

(A-6) 

are  solved  simultaneously  for  Up^. 

(A-7) 


A  "  Pof  ^t  ~  Pop  ^p  (A— 8) 

®  “  POf  ^t  2pOp  Sp  Vi  +  Pop  Cp  (A-9) 

C  -  -Pop  (Sp  V2i  +  Cp  Vl)  (A-10) 

If  the  projectile  and  target  are  similar  materials,  then  Eq.  (A-7)  be¬ 
comes  indetenalnant  and: 


Up  -  Vi/2  (A-ll) 

Once  Up^  Is  known,  Uj^  and  Pg^  ■  Pgp  may  be  computed  from  Eqs.  (A-2)  and 
(A-5),  respectively. 

The  problem  may  be  solved  graphically  via  plots  of  Pg  versus  Up  and 
Pg^  versus  Up^  as  shown  In  Figure  A-2.  The  graphical  solution  can  be^employed 
even  If  the  shock  velocity  Is  not  a  linear  function  of  the  particle  velocity. 
Either  Vi  or  Pg  can  be  given  and  the  graphical  procedure  will  supply  the 
unknowns . 


A-1 


I 


For  exADpla,  coniider  the  impact  of  iron  (or  mild  steel)  on  Comp'-BS  If 
the  Impact  velocity  is  2.02  km/sec.  From  the  x-marked  points  on  Figure  A-2: 

-  1.58  km/sec 

Un  •  0.44  km/eec 

Pp 

Pg  -  155  Kbars 

Convarsaly,  if  P,  •  155  Kbars  had  baan  glvan,  tha  particle  and  Impact 
valocltlas  could  also  ba  datermlnad. 

Tha  following  dtaenatonal  unit  relationships  proved  very  useful  in  the 
present  computations  and  literature  review: 


Ps 


Oo  ^S  ^p 


Gram  cm  cm 
cm^  y-sec  g-sac 


■  - Gr** -  «  Megabars 

cm  (g-sec)^ 

(A-12) 

P  (Kbars)  •  10^  P(Magabar8) 

■  10  P(Gigapascals) 

• 

-  10“®  P(Pascals) 

■  ^  (^) 

(A-13) 

P  P(Kbars)  -  10^  P(Mbar8) 

(A-14) 

1  bar  ••  14.504  Ibs/in.^ 

(A-15) 

•  0.98692  ataospharas 
sa  1  atmosphere 

Note  that  pressure  can  ba  expressed  in  units  of  energy  par  unit  volume 
(J/cm^) . 

Note  also  that  impact  shock  pressures  can  ba  enormous.  In  tha  above 
example , 

Pg  -  155  Kbars  -  (155)(14504)  PSI 
-  2,248,120  PSI 
2.25  Hillion  PSI 


A-2 


1.  Initial  shock  geometry  for  planar 


APPENDIX  B 


SHOCX  COMPRESSION  INFORMATION  FOR 
PROJECTILES  AND  EXPLOSIVES 


Moulard  [2]  employed  steel  projectiles  In  his  work  with'  Comp-B3.  Eldh 
and  coworkers  [8]  utilized  copper  projectiles  in  their  experiments  with  TNT- 
Slade  and  Dowey  [6]  used  both  brass  a. d  steel  projectiles  to  acquire  Impact 
detonation  data  for  both  Comp-B  and  Tetryl.  The  data  from  [6]  is  also  shown 
In  [7].  Steel  projectiles  were  employed  in  the  impact  detonation  experlMnts 
on  PBX  9404  reported  in  [1]- 

In  our  present  reexamination  of  these  data,  we  hava  used  shock  loaded 
Iron  data  for  the  steel  projectiles  which  were  assumed  to  be  mild  or  low  car¬ 
bon  steels.  For  the  brass  projectiles,  shocked  copper  data  was  employed.  The 
shock  compression  information  for  copper  and  iron  %rere  obtained  from  [10]  and 
[II]  which  are  standard  sources  for  such  Information.  The  data  used  are  s-iven 
in  Table  B-I .  See  Appendix  A  for  their  significance. 

Shock  compression  information  for  the  explosives  Comp-B3,  Comp-B,  TNT, 
Tatryl,  and  PBk-9404  is  given  in  Table  B-2.  The  data  sources  are  listed  also. 


Table  B-i.  Shocked  Projectile  Information. 


METAL 

Po 

S 

SOURCE 

Grams 

cm 

cm3 

U'sec 

Iron 

7.840 

0.3850 

1.580 

Refs.  10,  11 

Copper 

8.903 

0.3958 

1.497 

Refs.  10,  11 

Table  B-2.  Shocked  Explosive 

Information . 

EXPLOSIVE 

Oo 

Co 

SOURCE 

Grams 

cm 

cm3 

u-sec 

Comp-B3 

1.70 

0.303 

1.73 

Ref.  12 

Comp-B 

1.70 

0.295 

1.67* 

Ref.  12 

TNT  (Pressed) 

1.635 

0.208 

2.35 

Ref.  26 

(^*s  <  0.19 

0.020 

4.50  1 

Tetryl  (Porous) 

1.50  < 

1 

( 

Ref.  13** 

(Ug  >  0.19 

0.090 

2.67  ) 

P3X-9404 

1.842 

0.245 

2.48 

Refs.  14,  26 

\ 


*S  vas  modified  from  I. 58  to  1.67. 

**Ug  is  a  nonlinear  function  of  Up  for  pq  ■  1.50  grama/cm^.  Table  2  and 
Figure  3  of  [13]  vere  used  to  obtain  Ug,  Up  and  Pg.  This  corresponded  closely 
to  the  density  (1.54)  of  the  tetryl  employed  in  the  work  reported  in  [6].  The 
two  linear  relations  in  Table  B-2  vere  employed  in  the  present  analysis  to 
represent  the  nonlinear  variation. 
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As  mentioned  in  the  analysis  section,  the  critical  energy  criteria  w^s 
developed  by  the  authors  of  [4]  via  examination  of  experimental  data  reported 
in  [5].  These  data  were  for  explosive  detonation  via  thin  metallic  foil  im¬ 
pact  (Figure  12).  A  dependence  on  the  foil  thickness  (and,  hence,  shock  pulse 
duration)  was  noted  for  these  one-dimensional  impact  conditions.  Analysis  of 
these  data  revealed  that: 


^Pex 


P^s  ter 
(Po  Ua)ex 


-  Cl 

■  Energy/Area 
•  Force/Length 

-  Constant  (C-1) 


Where : 

■  Shock  Pressure 

■  Shock  velocity  of  the  explosive 

■  Particle  velocity  of  the  explosive  behind  the  shock 

Oq  ■  Unshocked  density  of  the  explosive 

t^r  ■  Constant  pressure  pulse  width,  2W/Ugp^^j 

w  «  Thickness  of  metallic  foil  or  place 

'^•proj  "  Shock  velocity  in  metallic  foil  inpactor 

C;  *  A  constant  or  practically  a  constant  for  many  explosives. 
It  has  the  units  of  energy  per  unit  area  or  force  per 
unit  length.  The  magnitude  of  Ci  Is  different  for  dif¬ 
ferent  explosives. 

Numerous  investigators  have  provided  comments  on  the  physical  signifi¬ 
cance  of  the  two  forms  (Pg  'p  t  and  Pg^t)  of  Eq.  (C-1).  Several  papers  pre¬ 
sented  at  the  Sixth  (1976),  and  Seventh  (1981),  International  Symposiums  on 
Detonation  ware  devoted  to  this  topic.  For  example,  see  [18]  through  [25] 
and  [40] . 

All  explosives  do  not  conform  to  Che  critical  energy  criteria  expressed 
by  Eq.  (C-1).  Roth,  in  [9],  discusses  the  implications  and  provides  examples 
for  explosives  which  follow  a  more  general  form  of  Eq.  (C-1)  which  la: 

p2gtcr  -  f(Ps)  (C-2) 

Whete  f(P|)  is  a  function  of  the  shock  pressur.. 
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In  this  connection,  it  should  be  noted  that  according  to  the  Cole  [36], 
the  fora  of  the  energy  flux  density  relation  (EqCC-l))  omits  the  contribution 
of  the  Internal  energy  (E-Eq)  and  kinetic  energy  (dp2/2),  and  considers  onl; 
the  pressure  work  (Pg/P).  If  the  Internal  and  kinetic  energies  are  accounted 
for,  theu  the  energy  flux  density  (energy  per  unit  area)  Is: 


(C-3) 


Since  In  the  present  case,  ell  quentltles  to  be  Integrated  are  constant 
from  t  •  0  to  t  ■  t(.r»  then: 


Kfd 


Pg  ^p  *^cr 


P*  Up  t(.f 


Thus : 
^s^  ^cr 


Po  Efd 


(C-4) 


(C-5) 


Which  Is  similar  In  form  to  Eq.CC-Z)  since  the  RHS  of  Eq.  (C-S)  Is  a 
function  of  Fg.  The  addition  of  the  Internal  energy  and  kinetic  energy 
appears  Important  enough  to  be  taken  into  account  In  some  cases  since: 


Varied  froa  1<13S  to  1.378  for  the  5  sets  of  data  considered  herein. 
That  Is,  the  Internal  and  kinetic  energies  contributed  froa  13  to  37  percent 
more  energy  than  accounted  for  by  Eq .  (C-1). 
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1 


This  would  appear  to  b«  a  cubjact  for  further  exploration  for  both  the 
critical  energy  and  critical  area  concepts,  particularly  for  those  explosives 
which  do  not  coiaply  with  Eq.  (C-1). 

The  fact  that  the  detonation  energy  criteria  for  oany  explosives  is  given 
by  Eq.  (C-l)  ioplies  that  the  shock  pressure  energy  (Pg/p)  is  the  doeinanr 
ignition  factor  for  these  explosives,  since  the  Internal  energy  (E  -  Eq)  and 
kinetic  energy  (Up^/2}  are  not  included.  Andersen  [41]  agrees  with  this 
interpretation  since  he  says: 

"Thus  the  critical  energy  nerely  represents  the  Hugoniot  energy  delivered 
during  the  shock  ignition  of  the  aaterial,  and  has  no  special  relevance  Co  Che 
IgniClon  process  ocher  chan  chat  a  particular  pressure  Is  required  to  cause 
self-sustained  Ignition  of  the  laaCerial  in  a  particular  tine." 

The  authors  of  [40]  provide  the  following  interpretation  of  Pg/p  work 
baaed  on  aulecul^r  dynamics  calculations: 

"It  appears  that  what  has  been  called  Pdy  work  is  really  not  a  piston 
action  but  a  suasation  of  Che  aicroacopic  kineaatic  processes." 

Also,  sea  Cole's  coivancs  ([36],  pp.  145,  233-234)  relative  to  shock  wave 
energy  dissipation  In  water. 

Noainal  values  for  the  constant  critical  energy  per  unit  area  (C],)  for 
Coatp-B3,  Coap-B,  TNT,  Tatty 1,  and  FBX  9404  (as  eaployed  in  the  present 
Investigation)  are. listed  la  Table  C-1.  The  sources  for  this  inforoatlon  are 
listed  also,  along  with  the  density. 

As  indicated  in  Tables  B-1,  B-2,  and  C-1,  the  shock  coapresslon  infor- 
aacion  necessary  to  analyte  the  iapact  detonation  data  cana  from  eultiple 
sources.  When  Information  from  two  or  more  sources  (even  for  Che  same 
designated  explosive)  must  be  eaployed,  there  is  still  considerable  risk  of  an 
"apples  and  oranges"  mix.  This  is  because  all  the  factors  which  affect  deto¬ 
nation  (coaposition,  aanufacturing  process,  projectile  aaterial,  grain  site, 
density,  etc.)  are  not  necessarily  the  seas  in  all  respects. 

For  exsBple,  Houlard  ([3],  Figure  2)  illustrates  differences  in  the  onc- 
dinansional  Pg-C  sensitivity  results  for  the  French  and  D.S.  Coap-B  explo¬ 
sives.  Apparently,  the  chemical  coapoeitions  are  slightly  different.  In  the 
present  analysis  of  Coap-B  data  [6],  a  value  of  Ps^pg^^c^  "  J/ca^  (Table 
C-1)  was  eaployed.  This  may  not  correspond  to  the  U.S.  Coap-B,  but  this  value 
froa  [15]  was  all  chat  was  available  when  the  iapact  data  ware  being  analysed. 

Moulard  ([3],  Figure  2)  also  notes  that  in  [2]  "Coap-B3  and  Cosip-B  1-D 
sensitivlcy  data  were  abusively  assuaed  identical."  If  he  is  referring  to  the 
French  InstiCute  of  Saint  Louis  (ISL)  Coap-B  and  Coap-B3,  there  does  not  seea 
to  be  ouch  difference.  We  were  able  Co  duplicate  Moulard's  [2]  results  by 
using  the  Information  listed  In  this  appendix  and  Appendix  B. 
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Although  It  could  b«  fortuitous,  the  information  listed  in  Tables  B-1, 
B'‘2,  and  C-1  gave  rather  consistent  and  reasonable  results  when  applied  to  the 
analysis  of  the  experimental  impact  detonation  data.  If  more  appropriate  in¬ 
formation  becomes  available  for  any  of  these  explosives,  the  basic  impact  data 
will  be  reexamined  to  ascertain  how  well  it  conforms  to  the  critical  area  con¬ 
cepts. 


TABLE  C-1.  Explosive  Critical  Energy. 


EXPLOSIVE 

PaUpg^tcr 

SOURCE 

Graaa 

Joulaa 

ca^ 

Cal 

cb2 

CoBp-B3 

1.73 

140.00 

33.4 

Raf.  2 

Coap-B 

1.73 

185.00 

44.2 

Raf.  15* 

TNT  (Prasaed) 

1.645 

142.36 

34.0 

Raf.  4 

TETBYL 

1.655 

42.00 

10.0 

Raf.  16 

PBX-9A0A 

1.842 

50.244 

12.00 

Raf.  4 

*Raporcad  in  a  ravlalon  Co 
caloria  •  4.187  Joulaa  Chan 

(15) 

cha 

.  Tha  priaary 
aoargy  ualta 

aourca  la  [3]. 
convaralon  factors 

Slnca  oaa 

ara: 

E(Joulas  ■ 

(4.187)  E(Calorlaa) 

- 

E(Calorlaa) 

-  (0.2388) 

E(Joulaa) 

C-5/(C-6  blank) 


Moulird's  critical  araa  concept  la  beat  explained  with  the  aid  of  Fig¬ 
ure  1  which  la  an  adaptation  of  Figure  2  of  [2]-  Figure  1  la  eaaentially  a 
more  detailed  depiction  of  the  shock  reflections  In  the  projectile  than  la 
shown  In  Figure  A-l  (Appendix  A)* 

Inspection  of  Figure  1  reveals  that  at  time,  t,  only  a  portion  of  the 
projectlle/exploalve  Interface  is  -fnlll  subjected  to  the  original  high  shock 
pressure  (Pg)  generated  at  impact  (t  -  o).  Only  that  portion  of  the  area 
enclosed  by  the  radius,  r,  has  been  loaded  with  the  initial  shock  pressure 
during  the  time,  t.  If  t  is  known  (t  <  R/Cp),  then: 

r  -  R  -  Cpt 

Moulard  defines  the  critical  Initiation  conditions  as  follows  (which  Is 
a  direct  quote  from  [2]): 

“At  a  given  shock  pressure,  the  shock-detonation  transition  occurs  only 
if  a  sufficient  area  of  the  explosive  target  is  loaded  at  high  pressure  for  a 
time  equal  or  greater  than  the  critical  shock  duration  t(p)  measured  in  plans 
shock-wave  experiments  at  the  same  shock  pressure.'* 

The  critical  shock  duration  time,  ter.  known  via  Eq.  (1)  If  either 
Pg  or  Vi  Is  known  (see  Appendix  A). 

^cr  ■  £L 

^^Pex 

The  critical  area  is  then: 

Acr  -  irR2cr 

Where: 

Rer  -  R  -  Cp  ter 

It  should  be  noted  that  the  critical  area  concept  is  independent  of  the 
constant  energy  (per  unit  area)  concept.  However,  t^r  «ist  be  known  from  some 
source  as  a  function  of  shock  pressure  or  some  known  variable.  For  exatsple, 
see  Eq.  (C-2)  in  Appendix  C. 

In  equation  D-4,  Cp  is  the  velocity  of  the  rarefaction  wave  from  the  edge 
of  the  projectile.  In  our  data  reduction ,  we  employed  Ug  ,  the  projectile 
shock  velocity,  for  Cp. 

The  above  procedure  to  determine  A^r  ^  adapted  for  different  projec¬ 
tile  cross-section  shapes.  Moulard  also  cleverly  employed  rectangular  and 
concentric  circular  cross-sections  to  illustrete  his  point.  We  were  able  to 
verify  his  calculations  of  A^r  for  these  shapes  (Tables  D-1  and  D-2).  These 
tables,  In  conjunction  with  Table  1  in  the  main  body  of  this  report,  illus¬ 
trate  the  procedure. 


(D-2) 
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The  A^j.  computations  for  the  tubular  or  annular  projectile  [2]  are  as 
follows  (Figure  D-2): 

R  -  Rl  -  Ugp  *  ter 

ri  “  r  ■  Ugp  *  ter 

Aer  ■  ir(Ri^  ~  ri^)  •  Are*  shocked  to  level  for  time,  ter* 

The  A^ji  computation  for  the  rectangular  bar  projectile  [2]  la  as  follows 
(Figure  D-3): 


/  L  -  Ugp  *  ter 
W  -  Ugp  *  ter 

Ll  -  L  -  2  (aL) 
Wi  -  W  -  2  (aW) 
Acr  -  Ll  *  Wi 


TABLE  [>-l.  Annular  Projectile  Critical  Area. 


ANNULAR 

PROJECTILE 

R 

r 

n 

Acr* 

Acr** 

A 

cm 

cm 

cm 

cm 

cm2 

cm2 

cm2 

A(9xl5) 

0.75 

0.45 

0.6817 

0.5183 

0.6160 

0.62 

1.131 

A(5xl5) 

0.73 

0.25 

0,6373 

Q.Zkll 

0.9883 

0.95 

1.571 

*See  Table  1 


‘♦Reference  2,  Figure  5 


TABLE  D-2.  Rectangular  Projectile  Critical  Area 


RECTANGULAR 

PROJECTILE 

L 

V 

Ll 

Wl 

Acr* 

Acr** 

A 

ca 

ca 

ca 

ca 

ca^ 

cb2 

ca2 

R(5xll) 

1.10 

0.50 

0.9916 

0.3916 

0.3883 

0.39 

0.53 

*S€e  Table  1 

**Refereoce  2, 

Fi gure  5 
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APPENDIX  E 


THE  JACOBS-aOSLUND  EQUATION 


An  caplrlcal  relation,  known  as  the  Jacobs  equation  or  the  Jacoba-Roalund 
equation  ([1]  and  [17])  is  widely  used  to  predict  the  critical  velocity  for 
projectile  induced  detonation  of  an  explosive.  This  equation,  in  a  practical 
fora  from  [1]  is: 

A  /  CT  \ 

®  @  d) 

Where: 

■  required  fragment  iispact  velocity  on  warhead  case  -  L/t 
A  •  explosive  sensitivity  coefficient  •  (L)^/^/t 
B  ■  fragment  shape  coefficient,  dimensionless 
C  ■  cover  place  protection  coefficient,  dimensionless 
T  ■  cover  plate  thickness  •  L 

0  ■  characteristic  dioMnaion  of  the  impactor  cross-'sectional  area  ■  L 

The  significance  of  the  factors  ®,  ®,  and  (J)  is  as  follows: 
factor  (l) 


VF 


V(.  for  impact  detonation  of  bare  explosive  with  flat  faced 
projectiles. 


Factor 


(1  -h  B)  ■  Correction  factor  for  projectile  shape  effects  (bare 
explosive) . 

Factor  (T) 


Correction  factor  for  cased,  covered,  or  protected 
explosive. 
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Nott  th«t  factors  and  are  both  concerned  with  the  bare  unpro¬ 
tected  explosive.  The  analysis^  this  report  Is  relevant  to  factor  1  for 
planer  lapact  with  flat  faced  projectiles. 

For  this  situation,  Eq.  (E-1)  reduces  to: 


(E-2) 


Cr 

v2j.d  ■  ■  Constant 

Note  that  Eqs.  (E-2)  and  (E-3)  are  equivalent  dimensionally  to: 


Energy 

Area 


l3v2 


(E-3) 


(E-A) 


and  p  1l  constant  If  the  projectile  oaterlal  reaalns  unchanged. 

Table  E-1  and  Figure  E-1  contain  the  computed  values  of  A  for  the  explo¬ 
sives  conaidared  In  the  present  study.  Figure  E-1  reveals  that 
indeed  constant  to  a  reiurkable  extent.  Table  E-2  lists  the  present  average 
results,  with  the  projectile  material.  Considering  the  factors  (see  Appendix 
C  discussion)  which  affect  sensitivity  for  a  given  explosive ,■ the  present 
values  appear  reasonable. 

The  eimplicity  of  Eqs.  (E-2)  and  (E-3)  can  hardly  be  improved  upon. 
However,  the  situation  is  more  complex  than  it  appears  to  be,  particularly  for 
noncircular  or  annular  cross-sections.  For  the  three  such  cases  considered 
hereir,  it  Is  worth  noting  that  computing  gave  results  consistent 

with  the  solid  circular  cross-sections  (Table  E-1).  D^q  is  the  diameter  of  an 
equivalent  circular  area. 

The  critical  area  concept  can  obviously  be  employed  to  check  the  "D” 
computations.  For  example,  the  projectile  cross-section  area,  A  must  be 
greater  than  A^^.  From  a  conservative  point  of  view,  the  smallest  dimension 
(W)  of  a  rectangular  projectile  is  the  most  iaiportant  one  from  Ug  t^^ 

consider itions .  The  example  of  the  R(S  X  11)  projectile  in  Appeolix  C 
illustrates  this.  Obviously,  W  must  be  greater  than  2Ug  t^.^.  or  A^f  is  zero 
and  detonation  Initiation  will  not  occur.  ^ 
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Table  E-1.  Detenainatlon  of  Explosive  Sensitivity  Coefficient. 


(Vc  VT  -  CONSTANT) 


EXPLOSIVE 

PROJECTILE 

D 

Vc 

VgV^ 

0m 

0m 

- 

m3/2 

u“8ec  y“aec 


C0MP-B3 

C(5) 

5.00 

2.020 

4.517 

C0MP-B3 

R(5xll) 

8.37* 

1.415 

4.093 

C0KP-B3 

A(9xl5) 

13.65* 

1.275 

4.710 

C0MP-B3 

C(10) 

10.00 

1.225 

3.874 

C0MP-B3 

A(5xl5) 

14.14* 

1.105 

4.155 

C0MP-B3 

C(15) 

15.00 

1.000 

3.873 

4.i04av 

COMP-B 

C(6.35) 

6.35 

1.320 

3.326 

COMP-B 

C(12.70) 

12.70 

0.900 

3.207 

COMP-B 

C(14,94) 

14.94 

0.850 

3.285 

TTITIv 

TNT 

C(6.00) 

6.00 

0.770 

1.886 

TNT 

C(7.50) 

7.50 

0.690 

1.890 

TNT 

C(IO.OO) 

10.00 

0.599 

1.894 

TNT 

C(12.00) 

12.00 

0.548 

1.898 

TNT 

C(15.00) 

15.00 

0.492 

1.906 

TS^Tav 

TETRYL 

C(3.175) 

3.175 

0.804 

1.433 

TETRYL 

C(7.620) 

7.620 

0.541 

1.493 

TETRYL 

C(12.700) 

12.700 

0.433 

1.543 

TETRYL 

C(14.940) 

14.940 

0.387 

1.496 

1 .491av 

PBX-9A04 

C(1.270) 

1.270 

2.20 

2.479 

PBX-940A 

C(1.778) 

1.778 

1.62 

2.160 

PBX-9404 

C(3.556) 

3.556 

1.13 

2.130 

PBX-9404 

C(7.620) 

7.620 

0.75 

2.070 

PBX-9404 

C(10.160) 

10.160 

0.65 

2.072 

PBX-9404 

C(19.050) 

19.050 

0.50 

2.182 

THTTav 

NOTE!  Vg  (■i/u“iec)  •  Vg  (k«/aec) 

*D  ■  square  root  of  cross-section  area  tlaes  2/V^ 
■  Equivi^rent  Diaaster,  D^q. 
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Table  E-2.  Explosive  Sensitivity  Coefficients  (A). 


PROJECTILE 

EXPLOSIVE 

MATERIAL 

A 

bb|3/2 

PRESENT 

RESULTS 

C0MP-B3* 

Steel 
(Ref.  2) 

4.204 

COMP-B 

Steel  and  Brass 
(Ref.  6) 

3.273 

TNT  (Pressed) 

Copper 
(Ref.  8) 

1.895 

TETRYL  (Porous) 

Steel  and  Brass 
(Ref.  6) 

1.491 

PBX-9A04 

Steel 
(Ref.  1) 

2.182** 

*Sce  cooaents  in  Appendix  C  vlth  respect  Co  this  explosive. 


**Re£ereace  1  deCenalned  A  to  be  2.05,  possibly  by  using  all  points  defining 
the  detonation  threshold  (Including  chose  Just  below  the  threshold). 
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E-5/(E-6  blank) 


Figure  E-1.  V^,  Y  D  versus  Vq  for  five  exploeiveB. 


APPENDIX  P 

DIAMETER  EFFECT  ON  DETONATION  StfOCE  VELOCITY 
FOR  CYLINDRICAL  EXPLOSIVE  SPECIMENS 


It  Is  well  known  that  there  is  a  size  effect  on  the  detonation  velocity 
achieved  lu  cylindrical  exploalve  samples  (for  example,  see  [27]-[33]).  That 
Is,  the  final  detonation  velocity,  D,  achieved  In  a  cylindrical  explosive 
tpaciaaa  (Ignlcsd  on  one  end)  Is  dependent  on  the  diameter  (d)  of  the  sample, 
normally,  there  Is  also  a  limiting  small  diameter  such  that  detonation  will 
not  propagate  (D  ■  0).  Equation  (1)  of  [29]  is  a  rather  general  form  of  the 
observed  relationships  between  D  and  d.  This  equation  Is: 


Where 


D«  ■  Detonation  shock  velocity  in  an  infinite  explosive  medium. 

D  ■  Detonation  shock  velocity  In  a  cylindrical  exploalve  sample  whose 

diameter  Is  d^^* 

^cx  "  Diameter  cf  explosive  specimen. 

R2L  ■  React icn  ts  Length,  or  some  characteristic  length  associated  with 
the  detonation 

The  authors  of  [29]  and  [30]  consider  the  reaction  zone  length  to  be: 

RZL  •=  (D  -  Up)  tr  (P-2) 

Where 


Up  «  Particle  velocity  In  the  reaction  zone  behind  the  detonation  shock 
front . 

tf  *  Reaction  tlma. 

Thus  Eq.  (P-1)  could  be  written  as: 

Wliere  f  still  denotes  a  general  functional  relationship  [31].  Both  theoreti¬ 
cal  and  experimental  results  for  many  explosives  indicate  that  bP  Is  directly 
proportional  to  l/d«x  explicit  functional  relationship  Is: 

-  k 

p«  L  J 


Where 


K  -  A  numerical  constant. 
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Physical  reasons  for  the  size  effect  on  the  detonation  velocity  are 
seated  by  Kistlakowsky  [27]  and  by  Ubbelohde  [28].  They  say  that  if  radial 
expansion  occurs  in  a  cylindrical  charge  behind  the  wave  front  before  the 
cheaical  reactions  have  been  completed,  the  cheticel  energy  of  this  unreacted 
explosive  la  uravailable  to  drive  the  wave.  The  amount  of  radial  expansion  is 
governed  by  the  diameter  of  the  charge;  the  aoialler  Che  diameter,  the  greater 
the  radial  expansion  and  chemical  energy  loss.  This  chemical  energy  loss 
lowers  the  detonation  shock  velocity  from  D„  to  D. 

Table  F-1  is  a  list  of  detonation  shock  information  for  Che  explosives 
analysed  in  this  report.  Sourr^s  for  this  information  are  given  also.  This 
information  is  from  one-dimensional  test  results. 

A  question  which  naturally  arises  is:  What  relation  does  the  phenomena, 
described  above  in  this  appendix,  have  Co  Che  topic  of  the  main  body  of  this 
report  where  the  explosive  samples  are  large  enough  (relative  to  the  projec¬ 
tile  size)  to  be  considered  an  infinite  tnediua? 

A  glance  at  Figure  1  (^hich  is  reasonably  scaled)  reveals  that  even 
though  the  explosive  tarpst  is  large  relative  to  the  projectile,  the  highly 
shocked  critical  region  where  Che  ignition  reaction  begins  is  confined  to  a 
volume  whose  dimensions  are  the  same  order  as  the  projectile  diameter.  Gener¬ 
ally,  as  shown  herein  (Tables  1  through  6),  these  dloeoalons  are  always  less 
Chan  the  projectile  diameter.  In  addition,  digital  simulations  which  model 
Che  molecular  dynasUcs  of  two-dimensional  iapacCed  solids  ([34]  and  [35]  show 
that: 

"Lateral  transfer  of  shock  energy  la  minimal,  even  when  Che  lattices  are 
initially  thermally  highly  excited." 

Thia  ia  a  direct  quote  from  [35]  and  aeaoa  that  most  of  the  shock  effects 
are  confined  to  the  material  Immadlatcly  ahead  of  the  finite  projectile. 
Numerous  figures  in  [34]  and  [35]  illustrate  this  also. 

Thus,  it  is  not  surprising  that  there  is  a  strong  projectile  ct^ss- 
seccion  size  effect  on  explosive  detonation  via  impact.  Furthermore,  certain 
aspects  of  finite  projectile  induced  detonation  may  be  aomawbat  similar  or 
analogous  to  certain  features  of  the  cylindrical  explosive  diameter  effect  as 
delineated  above  in  this  appendix. 

Sufficient  information  was  not  available  in  [1],  [2],  [6],  and  [8]  (or 
other  sources)  to  ascertain  if  Eqs.  (F-3)  or  (F-4)  were  applicbla  to  the  pro¬ 
jectile  Induced  detenation  data  specifically  considered  In  this  report. 
However,  somewhat  analogous  to  Eq.  (F-2),  an  ignition  zona  length,  IZL,  could 
be  computed  as  follows: 

IZL  -  (U,^^  -  ’Jp^^)  ter  (y-5) 

Where  L'g  ,  Up  and  t^r  hev«  the  same  definitions  employed  in  the  Min 
body  of  this  fiport**  See  also  Figure  1  which  shows  that  this  is  the  dis¬ 
tance  between  the  projectile  and  the  shock  front  when  t  -  t^r*  Ignition  phe- 
noMna  must  occur  within  this  length  of  shock  compressed  explosive  Mterial. 
Was  there  som  relation  between  the  IZL  and  the  projectile  cross-section 
dimensions? 
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To  answer  this  question,  the  ratio,  IZL/V^  computed  for  all  the  data 
coneidored  herein.  A  is  the  cross-section  area  of  the  projectile.  It  was 
found  empirically  that: 


IZL  .  (U»,^  -  Up.^)  ter 

yp: 


•  C5 

C5  is  a  constant  or  practically  a  constant,  over  a  wide  range  of  condi¬ 
tions  for  a  gtveu  explosive.  C5  was  different  for  each  of  the  five  explo¬ 
sives  considered  in  this  report.  See  Table  6  and  Figures  6  and  7  in  the  main 
body  of  this  report. 

The  ratio,  IZL/^,  is  soaewhat  analogous  to  the  ratio,  RZL/dgx*  00  the 
RflS  of  Eqa.  (F-1),  (F-3),-and  (F-4).  Equation  (F-6)  means  that  the  IZL  is 
directly  proportional  to  the  projectile  cross-section  area.  Note  that  this 
result  (Eq.  (F-6))  was  established  independently  of  Moulard's  Critical  Area 
Concept • 
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TABLE  F-1.  Explosive  Detonation  Characteristics 


EXPLOSIVE 

Po 

D. 

Up- 

"s- 

Graa 

ca 

ca 

eSAB 

c«3 

u~«ec 

ii~«ec 

C0MP-B3  1.70  0.790  0.280  376 


Ref.  12 

Ref.  12 

COMP-B 

1.70 

0.780 

Refs.  12,  33 

0.290 

Ref.  12 

385 

TNT  (Pressed) 

1.54 

0.690 

Ref.  31 

0.205 

Ref.  26 

219 

TETRYL  (Porous) 

1.50 

0.705 

Ref.  13 


PSX-9404 


1.84 


0.878 

Refa.  26^  33 


0.225 
Ref.  26 


360 


APPEliBIX  G 

MINIMUM  DETONATION  SHOCK  VAEIABLB  MAGNITUDES 
VIA  A  CRYSTAL  UTTICE  FRACTURE  CRITERION 


From  particle  dynamics  concepts,  Fitzgerald  (['^6],  Chapter  III)  derived 
an  expression  for  the  particle  velocity  necessary  to  produce  self-sustained 
lattice  disintegration.  He  called  this  velocity  the  phonon-fission  velocity, 

Vf 


Consequently,  it  seemed  feasible  that  lattice  breakup  could  initiate 
detonation  in  explosives.  If  so,  then  Vf  would  be  the  minlnum  particle  velo¬ 
city  (due  to  impact  or  ahock)  required  to  cause  detonation  in  one-dimensional 
large  samples.  The  corresponding  shock  velocity,  Ug^  and  shock  pressure,  Pgfi 
are  (Appendix  A) : 

Ugj  -  Co  +  S  *  Vf  (G-1) 

P*f  -  POe,  *  Ugf  *  Vf  (G-2) 

ELafercnce  [47]  documents  the  computation  of  Vs,  Ug^  and  Pg^  for  Comp-B3, 
Comp-B,  TKT,  PBZ-9404,  and  Tetryl.  Available  experimental  data  compared 
favorably  with  the  computed  results. 

Table  G-1  lists  Vf,  L'g^  and  Pg^  for  the  above  explosives  (47).  The 

Vf  magnitude  for  Tetryl  (po  ■  1.5  grams/cc)  le  considered  slightly  high  since 
certain  elastic  wave  velocities  (employed  to  compute  Vf)  were  for  oo  "  1*68 
grama/cc. 

Of  particular  interest  in  the  present  analysis  were  the  magnitudes  of 
Ug^  (lable  G-1)  which  would  be  employed  In  Eq.  (16)  of  Section  II. C  to  compute 
an  upper  bound  for  (2)tg>‘  Table  G-2  contains  the  nuawrical  results  for  both 
the  upper  and  lower  bounds  of  (E)t.f  os  per  Eq.  (16).  These  results  are 
plotted,  as  appropriate,  in  Flgurci^lZ,  13,  and  14.  See  the  discussion  in 
Section  II. C. 

For  comparison  purposes,  experimental  results  from  the  present 

analysis  were  acquired  via  Figures  G-1  through  G-S.  In  these  figures  the 
inverse  projectile  diameter  (1/dp)  is  plotted  versus  the  Initial  initiating 
pressure,  Ps.  True  one -dimensional  conditions  are  considered  to  exist  where 
1/dp  equals  zero. 

Note  that  for  these  explosives: 

(Pg  “  ^s^q)  *  Constant  (G-3) 

Pg^^  is  Che  extrapolated  poluc  where  1/dp  equals  zero,  or  one- 
dimensional  condiClona  axiac. 

Tha  computed  Pg^  and  experimental  Pa^Q  comparlaoQ  varies  from  fair  (TNI, 
Tatryl,  and  PBX-9404)  to  very  good  (Comp-B3  end  Coep-B). 

See  also  the  cheoreticel  end  exparlmentel  reeulca  comperieon  presented  in 
[47].  Ihla  illuetretes  Chet  the  megnltude  of  experlme  '■'.In  deCe  spread 
which  exlsce  for  ecriccly  one-dimenslooel  Cast  conditions  is  compsreble  Co  Che 
differences  In  Pg^  end  exhlbltsd  in  Figures  G-1  through  G-5. 


G-1 


EquaCioQ  G-3  Implies  that: 

(Pg  -  Pmin)  ■  Constant  (G-4) 

and  Eq.  (G-A)  aultipllad  by  Eq.  (6)  (Section  II. B)  yields: 

(Ps  ”  Pain)  ^^Sgj[  “  ^Pex^  *  Constant  (G—5) 

An  analyses  of  the  present  sets  of  projectile-explosive  impact  detonation 
data  has  conflraed  Eqs.  (G-4)  and  (G-5)  to  a  remarkable  extent*  See  Section 
II. B  in  the  main  text  of  this  report. 
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TABI.E  G-l.  NumerlcAl  Results  for  Vf,  Ugf  and  Pgf. 


MATERIAL 

Vf 

Usf 

"•f 

CB 

ca 

KBARS 

g-sec 

ij“sec 

CoBip-B3 

Po  -  1-70 

0.04170 

0.3751 

26.59 

Coap-B 

Po  ■  l-*70 

0.04077 

0.3631 

25.16 

TNT 

(Pressed) 

Po  -  1.635 

0.03739 

0.2959 

18.09 

TETRTL* 

(Porous) 

Po  m  1.50 

0.03543 

0.1794 

9.53 

PBX-9404 

Po  -  1.835 

0.04109 

0.3469 

26.16 

*Vf,  and  for  Tecryl  are  slightly  high  since  longitudlnel  (C]^)  and 
transverse  (C()  wave  velocity  data  for  "  1.68  graa/cc  were  eaployed  to 
coapute  Vf  [47 ] . 
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TABLE  G“2.  Computation  of  Upper  and  Lower  Bounds  on  (E)t|, 


EXPLOSIVE 

4 

TABLES 

1-5 

4 

TABLE  G-1 

4 

Poex 

D, 

TABLE  F-1 

4 

Oo  D_ 

^ex  • 

^Sf 

TABLE  G- 

Joulea^ 

J-CM 

CM 

Joules 

CM 

Joules 

CM^ 

X  106 

(u-Bec)2 

y-sec 

y-sec 

y-sec 

y-sec 

XBARS 

C0MP-B3 

41.5328 

244.31 

0.3751 

651.32 

0.7900 

309.253 

26.59 

COMP-B 

16.7690 

98.641 

0.3631 

271.66 

0.7800 

126.463 

25.16 

IN’!. 

0.2782 

1.7012 

0.2959 

5.75 

0.6900 

2.466 

18.086 

TETRYL 

0.2460 

1.640 

0.1793 

9.15 

0.705 

2.326 

9.53 

PBX-9404 

4.6708 

25.454 

0.3469 

73.38 

0.8780 

28.937 

26.16 

DESCRIPTION 


Inverse  prolectlle  diaweter  versus  impact  shock  pressure  for  Co«p-B3. 


TACLB  0-1,  P_  =  15.*  KBARS 


shock  pressure 


SYM  DESCRIPTION 


G-9/(G-10  blank) 


Figure  G-5.  Inverse  prolectlle  diameter  versua  Liapact  shock  preaaure  for  PBX-9404. 
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